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ABSTRACT 


Measuring of the masses of galactic supermassive black holes (SMBHs) is an 
important task, since they correlate with the host galaxy properties and play 
an important role in evolution of galaxies. Here we present a new method for 
measuring of SMBH masses using the polarization of the broad lines emitted from 
active galactic nuclei (AGNs). We performed spectropolarometric observations 
of 9 AGNs and find that this method gives measured masses which are in a good 
agreement with reverberation measurements. An advantage of this method is 
that it can be used to measure the masses of SMBHs in a consistent way at 
different cosmological epochs. 

Subject headings: (galaxies: nuclei quasars: emission lines - line: profiles - 
polarization 
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Introduction 


It is widely accepted that all massive galaxies host supermassive central black holes 
(SMBHs) ranging in mass from less than a million solar masses to many billions of solar 
masses. The masses of SMBHs are known to correlate with the bu lge properties of their 


host galaxies ([ Heck man & Ka uffm ann 


2011 


Kormendv & Ho 


this is important for understanding the evolution of galaxies (see 


20131). and understan din; 


Heckman fe Best 


20R 


Therefore, measuring of the SMBH masses at different cosmological epochs is an important 
task in astrophysics today. There are several methods (direct and indir ect) w hich have 


been used to measure black hole masses in the center of galaxies (see 


Peterson 


2014 1. where 


direct methods, especially reverberation, can be used for measuring black hole masses of 
low redshift quasars. However, measuring of SMBH masses from a single epoch observed 


spectrum of quasars is still in a developmental stage ([Peterson 


2014). 


Here we show how the polarization of the broad lines emitted from active galactic 
nuclei (AGNs) can be used to measure the masses of supermassive black holes. We 


tested the method (see 


Afanasiev et al 


2014) for a sample of nine AGNs which have been 


observed with 6-m telescope of Special Astrophysical Observatory of Russian Academy of 
Science (SAO RAS). In difference with reverberation, the method only needs one epoch of 
observation of an AGN. 

The paper is organized as following: in §2 we give method and in §3 we outline the 
basic of observations and obtained results. 


2. Method 

We make use of the fact that polarized light in broad AGN lines, can, in some cases, 
be interpreted as being caused by scattering in the inner part of a dusty torus (so called 
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equatorial polarization, for more details see lSmith et a,1 .1120051: lAfanasiev et al.l 120141) . The 
broad lines are emitted from the broad line region (BLR) that is close to the SM BH, 


consequently one can expect near Keplerian motion of the emission gas in the BLR (jGaskell 


20091 1. A rotating, Keplerian line-emitting BLR surrounded by a co-planar scattering region 
produces a polarized broad line, where the position angle (PA) of polarization avera ged 


across line profile is aligned with the projected BLR rotation axis flSmith et ah 


200511 . and 


the stratified velocity held of Keplerian-like motion in the BLR produces a characteristic 
change of the PA across the line profile (see Fig. [U upper panel). This gives nearly 


equal but opposite m aximum in the blue and red wings of the line (Smith et ab 


Afanasiev et al. 


2014). 


2005 


Let us consider a simple model as shown in Fig. |Tj (panel down), considering 
approximation of a single scattering element from the torus. The polarization angle 
across t he line profile, A<p(A) = <p^(A) — <^c(^)El, depends on the velocity held in the BLR 


(Smith et ah 


200511 . For Keplerian motion the velocity depends on the distance of the 


_1/2 

emitting gas from the SMBff as Vj ~ ' and, it also depends on the polarization angle 

as Vj ~ tan(A<pj) -1 / 2 (see Fig{H lower panel). The relationship 


Detween velocities and 


polarization angles across the line prohle is thus flAfanasiev et al. 


2014): 


log(—) = a - 0.5 • log(tan(A<pj)), 
c 


( 1 ) 


where c is the velocity of light. The constant a directly depends on the black hole mass as 

,GM bh cos 2 (9) 


a = 0.5 log(- 


c 2 R, 


( 2 ) 


where G is the gravitational constant, R sc is the distance of the scattering region from 
the central black hole and 6 is the angle between the disc and the plane of the equatorial 
scattering region (see Fig. [[]). Since the BLR is expected to be nearly co-planar with 


1 where <Pl(A) is the PA in the line and tpc/A) is the PA in the continuum 
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the torus, one can take 9 ~ 0 as a good approximation. The effect of a wide (or a non 
co-planar) torus (e.g. 6 ~ 10 — 20°) can give an error in black hole mass (see Eq. (2) 
M BH (6 = 0 )/M bh (9 % 0) = cos 2 (6>)) measurement around 5-10%. There also can be 
partial obscuration of the BLR, but it is connected with the orientation of the system 
with respect to the line of sight. In principle one cannot expect a high contribution of the 
equatorial scattering of the line light in the edge- or face on orientation of the system, i.e. 
the equatorial scattering can be expected in the systems with inclination between 20-70 
degrees. On the other hand, the relation between velocities and A ip does not depend on 
the inclination of the system, since the Keplerian disc emits nearly edge-on orientated line 
light to the scattering region and relation between log(R) and logftan(A gA) is close to 0.5 


(as it shows a preliminary simulation by STOKES code 


Goosmann 


2014). 


Note here that in the case where full scattering ring is taken into account (not a single 
scattering region as we assumed here), scattered emission will occur at various polarization 
angles and there will not be one to one relation between R and <p. However, in the case of 
inclined system (torus + Keplerian disc), the dominant scattered light is coming from one 


side of the torus, while scattered li, 


ght for opposite side is blocked, and the single scattering 


region approximation can be used flGoosnmnnl 2014). More detailed discussion of the model, 


numerical simulations and exploration of the BLR and scattering region parame ters us ing 
polarization in broad lines will be given in an extensive forthcoming paper (Afanasiev et al. 
20151). 


To measure the black hole mass using Eq. (2) it is 


expected to be the inner radius of the torus flSmith et ah 


n eces sary to estimate R sc . This is 


2005). There are two possible ways 


to find R sc . The first is to use the empirical re. 


luminosity of the AGN fjKishimoto et al 


2011 


ationship between th e inner torus radius and 


estimate it directly from infrared observations (Kishimoto et al 


Koshida et ah 


201 41), and the second is to 


2011 


Koshida et al 


2014|) 
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3. Observations and Results 


To test the method here we have selected a sample of nine AGN ( see Table CD), for 


whic 


1 estimates of t 


2011 


Koshida et ah 


re inner radius of their tori are given in literature (Kishimoto et al 


2 0140 . We have obtained spectropolarimetry of the AGN with 6-m 


telescope of SAO RAS using a modified version of the SCORPIO spectrograph (for more 


details, see 


Afanasiev fe Moiseev 


201ll) . The reduction of observed d ata and corrections 


for the interstellar polarization is described in 


Afanasiev et ah 


(20141) and here will not be 


repeated. The observed AGN and their basic data are given in Table [D 


The observed shapes of Atp and — for two of the objects, 3C273 and NGC4051, are 
shown in Fig. [21 It can be seen in Fig. [2] that the polarization angle shapes indicate 
Keplerian-like motion (Fig. [21 upper panel) and that Eq. (1) fits the observed points 
relatively well except for the low-velocity part (see Fig. [21 lower panel), where emitting gas 
is close to the torus and possibly Keplerian motion is not dominant. Therefore, we did not 
considered points with Vj ~ 0, since they are close to the inner radius of the torus and this 
simple scattering geometry does not work, then we take Knm for A ip < 45° 


In Table [T| we give estimat es of R sr with corresponding literature refe rences (columns 


4th and 5th, respectively - see 


Kishimoto et al. 


2011 


Koshida et ah 


20141 ). estimates for 


coefficient a (see Eq. (2)), and our black hole masses (columns 6th and 7th, respectively). 


We have also compiled reverberation- mapping estimates o f the black 


lole masses from the 


literature (column 8 in Table 1 - see 


2010 


Grier et al. 


Peterson et al 


2004; 

Benz et al. 

2007; 

Dennv et al. 


20121) and we compare these with our spectropolarimetric estimates (see 


Fig. ED. As it can be seen in Fig. El there is a good agreement between our measurements 
and reverberation ones. This supports the validity of our proposed method and the good 
correlation between line polarization and reverberation mapping masses strongly supports 
our assumption that the scattering region is at or near the inner radius of the torus. 
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Our method offers a number of advantages over traditional reverberation mapping (see 


Peterson 


2 014 ). The first and most fundamental advantage that is that our estimates of 


SMBH masses do not depend on the BLR inclination and geometry, i.e. the additional 
effects to the rotation of the BLR, as outflows/inflows, which may strongly affect the broad 
line profiles (especially widths, which have been used in the reverberation method), do 
not change the PA shape across the line profile that is used in this method. These effects 
can be seen only as velocity shift of the PA center (see Fig. [2j) . A second advantage is 
that our method only needs one epoch of observation whilst reverberation mapping is very 
telescope-time intensive. A third advantage is that while in the reverberation method it has 
to be assumed a priori that the BLR is virialized, in using the spectropolarimetric method 
one can test the assumption Keplerian motion (i.e., virialization) using the relationship 
between A <p and velocities across the broad line profile (see Eq. (2) and Fig. [2]). Finally, 
another significant advantage of our method is that in principle it can be applied to lines 
from different spectral ranges going all the way from the near infrared and optical (such 
as the Balmer lines used as in this paper) to the UV part (e.g., Lycq CIII], C IV and 
Mg II) so long as the variation in PA of the polarization can be measured. This thus 
allows measurement of black hole masses in a consistent way at different cosmological 
epochs. Doing this in practice requires the development of high quality spectopolarimetric 
instrumentation. 


On the other hand, there are some problems with application of the method to 
some AGNs: i) The inner torus radius measurements is needed for the black hole mass 
estim ates - that can be find by using the reverberation method in the infrared (2-3 /.irn, see 


Kishimoto et ah 


2 0111 ) for low-redshift AGNs, and using the calibration b etwe en the inner 


torus radius and the UV radiation for high-redshift quasars (see e.g. 


Barvainis 


1987); ii) the 


method can be used only for a rotating BLR disk, and in the case of the BLR with dominant 
radial component (without significant rotation), the method cannot be used. It may be a 













problem with high ionized lines, as C IV and CIII], however, there is indications that also 


in high ionized lines a Keplerian motion component is present (see e.g. 


Clavel et al 


1987). 


In near future we are going to perform new spectopolarimetric observations of high redshift 
quasars in the UV lines. 


Spectropolarimctric observations were obtained with the 6-m BTA telescope of the 
Special Astrophysical Observatory Academy of Sciences, operating with the financial 
support of the Ministry of Education and Science of Russian Federation (state contracts 
no. 16.552.11.7028, 16.518.11.7073). The authors also express appreciation to the Large 
Telescope Program Committee of the RAS for the possibility of implementing the program 
of Spectropolarimetric observations at the BTA. This work was supported by the Russian 
Foundation for Basic Research (grant N12-02-00857) and the Ministry of Education, Science 
and Technological Development (Republic of Serbia) through grant 176001 (“Astrophysical 
Spectroscopy of Extragalactic Objects”). L.C. Popovic thanks the Alexander von Humboldt 
foundation for supporting his research in general and COST Action MP1104 “Polarization 
as a tool to study the Solar System and beyond” for supporting his research in this field. 
We thank Martin Gaskell for useful discussion and Rene Goosmann for giving some 
explanation with unpublished model simulations. We thank to an anonymous referee for 
very useful comments. 
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Line-of-sign Scattering region 



Fig. 1.— Our assumed scattering geometry in indicated in the lower part of the figure and 
the observed dependence of polarization angle fA( p) vs. velocities in the Ha line profile of 


Mrk 6 is shown in the upper part of the figure (see 


Afanasiev et ah 


2014). 
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3C273, Rsc = 953 ligth days 



Log(V/c) = a + 0.5 Log(tamp) 


NGC4051, Rsc = 38 ligth days 



Log(V/c) = a + 0.5 Log(tan^) 




Fig. 2.— The observed shape of (Ay>) (upper panels) and velocities (lower panels) across 
the H line profile and for 3C273 (left hand) and NGC 4051 (right-hand panels) for the AGN 
in our sample with the highest and lowest masses, respectively. The full circles are from the 
blue and open are from the red line part. The gray curves in upper panel are expected Ay> 
for estimated black hole masses. 
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Fig. 3.— Comparison between black hole masses obtained by onr spectropolarimetric method 
with reverberation estimates in the literature (see Table [T]) . 
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Table 1: The list of observed AGN with basic data and measured black hole masses using the 
polarization in broad line (Column 7) compared with ones me asur ed by the reverberation 


method (Column 

8). T 

2- 

Koshida et al. 

2014) 


from 3 


Grier et al 


’lie estimates for Rsc- are taken from 1 -lKishimoto et al .1 (2011) and 
(I2014J1 . Th e reverberation meas ure ments of black hole masse s are taken 


(2007) 


(2012); 4-Peterson et al 


(2004); 5- 


Denny et al. 


(120101 ) and 6- 


Benz et al. 


Object 

Type 

z 

Rsc [pc] 

Ref. 

—a 

log (M pol )[Mq\ 

log (M rev )[Mq\ 

Ref. 

Mkn6 

1.5 

0.0188 

0.185 

1 

2.19T0.21 

8.18T0.42 

8.13T0.04 

3 

3C273 

1.0 

0.1583 

0.809 

1 

2.19T0.13 

8.85T0.27 

8.95T0.09 

4 

Aknl20 

1.0 

0.0323 

0.380 

1,2 

2.44T0.18 

8.02T0.36 

8.18T0.06 

4 

NGC4051 

1.0 

0.0024 

0.032 

1,2 

2.90T0.09 

6.23T0.18 

6.24T0.13 

4 

NGC4151 

1.5 

0.0033 

0.037 

1,2 

2.34T0.13 

7.21T0.27 

7.12T0.05 

4 

Mkn335 

1.2 

0.0258 

0.119 

2 

2.48T0.07 

7.44T0.14 

7.40T0.05 

3 

NGC3227 

1.5 

0.0039 

0.021 

1,2 

2.32T0.21 

7.01T0.42 

6.88T0.10 

5 

NGC5548 

1.5 

0.0172 

0.096 

2 

2.30T0.27 

7.70T0.33 

7.82T0.02 

6 

Mkn817 

1.5 

0.0315 

0.151 

2 

2.42T0.11 

7.67T0.21 

7.64T0.11 

5 






























